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Abstract 

The technique of decomposing Feynman diagrams at the one loop 
level into elementary integrals is generalized to the imaginary time 
Matsubara formalism. The three lowest integrals, containing one, two 
and three fermion lines, are provided in a form that separates out the 
real and imaginary parts of these complex functions, according to the 
input arguments, in a fashion that is suitable for numerical evaluation. 
The forms given can be evaluated for arbitrary values of temperature, 
particle mass, particle momenta and chemical potential. 
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I. INTRODUCTION 



The imaginary time Matsubara formalism for Green functions Q] is a venerated 
method of dealing with problems of finite temperature. This is so because of the 
fact that Wick's theorem can easily be shown to hold for the Matsubara operators. 
As such, this method of dealing with finite temperature systems has found wide 
application, primarily in condensed matter physics ||, but also recently in problems 
relevant for high energy physics. In particular, it has been recently used in studies 
of effective QCD Lagrangians (see, e. g. the reviews || and [Q and references 
cited therein), through which an understanding of hot and dense matter is sought 
[|§-£§. However, to date most such studies have often been constrained by specific 
parameter choices to enable tractability - often degenerate particle masses, zero 
chemical potential or special kinematics are chosen. Furthermore, often only the 
principal values of these integrals are considered, with the complex nature being 
completely ignored. An extension of those calculations to physically relevant cases 
with arbitrary parameters and generalized kinematics, as may be appropriate for 
example, in calculating cross sections in an SU(2) or SU(3) effective model of QCD 
10-01, becomes technically extremely complicated. It is thus the purpose of this 
paper to address the technical problems associated with such integrals, and present 
them in a form suitable for the numerical evaluation of both their real and imaginary 
components. This work has come about through a study of hadronization processes 
within the Nambu-Jona-Lasinio model ||, and is intended to provide the technical 
knowledge required for such functions for researchers who are active in this field. 
The associated computer code conformal to the analysis presented here is freely 
available?]. 

For calculations at zero temperature and density, it has been shown [|l(J, that 
all occurring transition amplitudes at the one loop level can be decomposed into a 
number of elementary integrals, which in turn can be classified by the number of 
particle lines they contain. Since this decomposition technique can be generalized 
to the Matsubara formalism, the aim of this paper is to provide suitable forms for 



Programs are available by anonymous ftp to Trick. MPI-HD.MPG.DE, the location of the 
archive file is /pub/loopies-1 . . shar . Z. 
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the elementary integrals containing one, two and three fermion lines, which we term 
A, B and Co, respectively, which can be used for numerical evaluation. (The cal- 
culation can be easily taken over to boson loops on replacing the Fermi distribution 
by the Bose distribution.) We analyze these integrals, explicitly displaying their 
complex nature that is a function of the input parameters. The parameters that 
may then be specified can be of the most general form, i. e. arbitrary tempera- 
tures, particle masses, chemical potentials and kinematics. As such, the associated 
computer code can be treated as a library-like 'black box'. 

We note at this point that the calculation of the one loop integral A is straight- 
forward. Its evaluation here serves merely to remind the reader of the usage of the 
Matsubara formalism and it also renders this work complete. The evaluation of B , 
on the other hand, serves to demonstrate the complexities involved in a consistent 
evaluation of this function. Finally we deal with the three line integral C , which 
then is highly nontrivial and complicated. 

It is shown in this paper that each of integrals B and Co (and thus implicitly all 
higher one loop integrals) can be written as a sum over terms that are constructed 
from a single generic function that is taken at different values of its arguments. 
The analytical structure of these generic functions is connected with the number of 
nonvanishing external momenta. We are thus required to distinguish between several 
cases in evaluating the generic functions when constructing a physical quantity. 

As is usually the case in relativistic quantum mechanics, the one and two fermion 
line integrals A and B are divergent, which means that they have to be regularized. 
We do this using a three momentum cutoff A, which means that we integrate over 
all momenta with \p\ < A. The three fermion line integral Co is however convergent. 
We nevertheless perform our analysis of this function including a three momentum 
cutoff. The reasons for doing this are twofold: (i) the extension to the case A — > oo 
is simple to obtain and provides a check of our calculations, and (ii) for models such 
as that of Nambu and Jona-Lasinio extended to QCD, it is regarded as internally 
consistent within the model to restrict all quark momenta in all loops. 

Having dealt with the analytical complexities of the loop integrals, we demon- 
strate in the last section the usage of the integral decomposition in a simple physical 
example, viz. that of calculating the kaon masses in the 577(3) Nambu- Jona-Lasinio 
model. This involves the integrals A and B . Use of C may be found, for example, 
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in Ref. 0. 

This paper is organized as follows: in Sec. we define the functions A, B 
and Co- These functions are brought into a form suitable for numerical integration 
in Sees. [TTI|, [TV] and M, respectively. At the end of Sec. [TV], we also give some 
technical hints for the evaluation of our results. In Sec. [VI], we give an example for 
the application of our calculations. Further examples can also be found in Ref. |§. 
In Sec. |VII| we summarize and conclude. Some special formulae are proven in an 
appendix. 



II. DEFINITION OF A, B AND C 

We define the functions A, B and Cq to be the analytical continuation of the 
one, two and three fermion line integrals, 

167T 2 f d^p 1 

A(m, „, A) = — Eexp(i^) J ^ + ^ _ ^ , (2.1) 

B (k, mi, ji h m 2 ,/i 2 ,w m , /3,A) = (2.2) 
167r 2 ^ , . s f d 3 p 

and 

C {k,q,6^, mi,fj,i, m 2 ,/i 2 ,^ m , m 3 ,fi 3 ,ia h (3,A) = (2.3) 



2^ exp K %u nm j (2?r)3 + _ _ + _ 



' 



167T 2 ^ /• rf 3 p 1 



(3 ^ ^ nu J (2ir)*((iu; n + fii) 2 -El)((iu; n -iv m + fi 2 ) 2 -El) 

1 



x- 



((iu;„-ia ; +/i 3 ) 2 -£ 2 ) 
in which the complex frequencies %v m and ia\ are to be continued to their values 
on the real axis after the Matsubara summation on n is carried out. It is to be 
understood that the limit rj — > is to be taken after the Matsubara summation. In 
the above expressions, we have introduced the abbreviations 



E = \Jp 2 + m 2 Ei = yp 2 + m 2 

(2.4) 



E 2 = V (P _ k ) 2 + m 2 E 3 = v (p - q) 2 + ™ 
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where the symbols and /x, denote the masses and chemical potentials of the 
fermions running on the individual lines of the diagram. The rrii are to be understood 
to include a negative infinitesimal imaginary part, mi — *> rrii — ie, e > 0. In Eqs. ( |2.1|) 
to ( |2.3| ), (3 is the inverse temperature, (3 = 1/T. The three momentum cutoff is 
denoted by A, and all integrals are considered for \p\ < A. 

The complex four momenta (iv m \ k) and (ia^ q), are the four momenta that enter 
the loop. The frequencies iv m and iai are bosonic in nature, 

2mir 21tc 

while the frequencies iuo n are fermionic, 

(2n + l)n 

u n = . (2.6) 

After analytical continuation, %v m and iai become the zeroth components associated 
with the three momenta k and q, i. e. iv m — > k , iai —> Qo- O ne may note at this 
point that, due to rotational invariance, the functions B Q and Cq do not fully depend 
on both k and q; B depends only on k = \k\, while Co depends only on k, q = \q\ 
and the enclosed angle 5^ as is implied by the arguments specified on the left hand 
side of Eqs. (|2.2|) and 



III. THE CALCULATION OF A 

The calculation of A is simple and straightforward. After evaluating the Mat- 
subara sum by contour integration in the usual fashion [Q], one has 

d 3 p 



A(m,fi, @,A) = 16^ 



(27T) 



4 [ AE dE VE 2 - m 2 (f(E -fi)- f(-E - //)) 

Jm 



where we have introduced the Fermi distribution function 

1 



exp([3x) + 1 



and the energy cutoff 



(3-1) 



(3.2) 



a e = Va 2 T 



(3.3) 
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From Eq. (|3.1|) , it follows that A(m,fi, 13, A) = A(m, —fi, /3, A). This integral can 
be easily performed by numerical integration. At T = 0, it has a closed analytical 
form, 



A(m, fi, oo, A) = 29(A| - ;U 2 ) ( m 2 log ^^.f + kk' - AA E ) (3.4) 



where k = max(m, and = k 2 — m 2 . 

In Fig. [TJ, we show A/ A 2 as a function of m/A, for the specific parameter choice 
/i = 0, (3 -> oo and A = 602.3MeV. At m = 0, we obtain the limit A(0, 0, oo, A) = 
— 2A 2 . For larger values of m, A increases continously to zero. 

IV. THE CALCULATION OF B 

A. Matsubara Summation 

The two line integral B (k, rax, fix, Tn 2 , [i 2 , ^m, (3, A) can be analyzed in a similar 
fashion to A(m,/j,, (3, A). The denominator of Eq. (|2.2j ) has the four zeros 

iu n = -/ii ± E 1 iu n = iv m - Hi ± E 2 , (4.1) 

from which we obtain, after the summation over n, the terms 

B (k, mi,/ii, m 2 ,p, 2 ,iv m , (3, A) 



16tt 2 



d 3 p 
(27p 



+ 



2Si (-//i + E x - iv m + /i 2 ) 2 - E 2 

H-Ex-fr) 1 

2Ex (-Mi -Ex- %v m + fi 2 ) 2 - E 2 

f(E 2 -» 2 ) 1 



2E 2 (iu m - /i 2 + E 2 + nx) 2 - Ef 
f{-E 2 -» 2 ) 1 



(4.2) 



2E 2 (iu m -ii 2 -E 2 + ii 1 ) 2 -Ef 

To eliminate the angular dependence in the arguments of the Fermi functions, we 
carry out the substitution 

p — > k — p (4.3) 
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in the third and fourth terms of Eq. ( |4.2|) and obtain 

B Q (k, mi,/ii, m 2 ,ii2,ivm, P,A) 



in 2 



d 3 p 
(2^)3 



J '(Ex -//i) 1 



+ 



Ei A 2 - 2XEi + 2pA; - k 2 + mj - m\ 

f(-Ei - in) 1 

Ei A 2 + 2A£i + 2pk -k 2 + ml-m 2 2 

f(E 2 -fi 2 ) 1 



E 2 A 2 + 2\E 2 + 2pk -k 2 -m\ + m\ 

_ f(-E 2 - /i 2 ) 1 

E 2 A 2 - 2XE 2 + 2pk -k 2 -ml + m 2 2 _ 
where the meaning of Ei and E 2 has changed to become 



(4.4) 



Ei = ^p 2 + m\ E 2 = ^p 2 + m 2 2 . (4.5) 

We have also introduced 

A = %v m + Hi - \i % . (4.6) 

Note that A will be real after the analytical continuation. One should further note 
that the substitution Eq. ( |4.3Q has been performed under the assumption that A 3> 
k. Thus no corresponding shift of variables in the integration limits has been carried 
out at this point. 

Without loss of generality, we may choose the coordinate system in such a way 
that k points in the z direction. Then the integration over the angle <fi yields a factor 
27T. It is now useful to consider each term occurring in Eq. ( [4.4| ) separately. These 
correspond to the individual poles in the original expression and have a common 
generic structure. One may therefore decompose Bq into the four terms 

B (k, mi, jji, m 2 ,jj 2 ,iu m , p, A) = S+(-A, k, mi,m 2 ,jji, (3, A) 

- B (\,k, mi, m 2 ,//i, (3, A) 
+ B^(X,k, m 2 ,mi,fi 2 , (3, A) 

- Bq(-X, k, m 2 ,mi,n 2 , (3, A) (4.7) 
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in terms of the generic function 
B (X,k, m,m',[i, (3,A) = 

Me r+l 1 

2 dEpf(±E-u) dx— (4.8) 

J m yJK J-i X 2 + 2XE + 2pkx - k 2 + m 2 - m' 2 V ; 

with p = E 2 — m 2 . Numerical computation of B thus relies purely on the eval- 
uation of B^. The analytical structure of this integral differs for k = and k > 0, 
so that these two cases must be considered separately. We examine these in the 
following two subsections. We may also note that the symmetry relation 

B (k, mi,/ii, m 2 ,/i 2 ,-A; , (3, A) = B (k, m 2 ,/i 2 , m 1 ,fx 1 ,k , f3,A) , (4.9) 

can be inferred from Eq. ( |4 . 4j ) . 



B. Calculation of B% for k = 



For k = 0, the integral over x is trivial, and we obtain 



B °- A L dE X 2 + 2XE + m 2 -m> 2 " ^ 10) 

(Here and in the following, we drop the arguments of Bq for convenience.) The 
analytical continuation of the Matsubara frequency iv m to real values leads to poles 
in the integrand at E = Eq, where 

X 2 + m 2 - m' 2 

E = - x , (4.11) 

if m < Eq < Ke- To calculate the integral in this case, we recall that the particle 
masses are complex, m 2 — > m 2 — ie, m' 2 — > m' 2 — ie, and apply the formula 

1 1 

lim = V- + i7r5(x) , (4.12) 

£->o x - ie x 

where V denotes the Cauchy principal value. From this, we find 



e-o Jm X 2 + 2XE + m 2 - m' 2 - zesgn(A) 



e->U A Jm hi H — ^ ie 



12 



X 2 + 2XE + m 2 -m 
+ i—p f{±E -f Jl )G{{A E -E ){E -m)) , (4.13) 



where Po = y Eq — rn 2 and the © function ensures that the imaginary part appears 
only if Eq lies within the integration interval [m, A E ]. 

Figure shows B as a function of fc /A for the parameter set /c = 0, /ii = /i 2 = 0, 
m i — m 2 — m — 367.7MeV, (3 — > oo, A = 602.3MeV. One recognizes that at small 
fco, -Bo is a real function, whereas for ko > 2m an imaginary part emerges. At 
k > 2A E , the imaginary part vanishes due to the presence of the cutoff. 



C. Calculation of Bq for k > 

For k > 0, the generic function 5^ has the form 
Bn(\,k, m,m r ,fi, /?,A) = 21im dE pf(±E - //) 

r+1 1 

x / dx , (4-14) 

J-i A 2 + 2\E + 2pkx — k 2 + m 2 — m' 2 — iesgn(A) ' 

where we have already indicated that the infinitesimal imaginary part in the inte- 
grand is determined by the sign of A. In the limit e — > 0, the integrand contains a 
singularity at 

1 



if 



x = — (A 2 + 2XE - k 2 + m 2 -m' 2 ) (4.15) 

2pk 



X 2 + 2XE - k 2 + m 2 -m 2 \ <2pk . (4.16) 



On performing the x integration and applying Eq. ( |4.12 ), one obtains 



^ = H \ t E dE f ( ±E - /*) r dx i x^2XE-W- m - ■ 7TT ( 4 - 17 ) 

e^o k Jm J-i x + + 2 , + — xesgn(A) 



f AE dEf(±E-/i) log 



l r A E 
k 



2pk 

\2 ( m j>\2 ™/2 



(A + Ef - (p - A;) 2 - m ,: 



(A + £) 2 - (p + k) 2 - m 



12 



+ ^7rsgn(A) f A \ E Q( 2pk _ lX 2 + 2XE _ k 2 + m 2_ m ,2 l)f{±E _ fi) 

fc Jm 
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The integral for the imaginary part can be done analytically: 

dE eM±PE- M + 1 = E T \ '° S(1 + eM±0E ~ M) ' 



(4.18) 



The integrand for the real part in Eq. ( |4.17|) contains logarithmic poles at the 
energies E = Eu% which are defined by 



\X A + 2XE 1/2 - k' A + m z - m a \ = 2k J E 



J l/2 



(4.19) 



These energies also constitute the limits of integration for the imaginary part. Squar- 
ing Eq. (|4.19|) leads to a quadratic equation for Ex/2, which has the solutions 



E 



1/2 



A(A 2 - k 2 + m 2 - m' 2 ) , J(A 2 - k 2 + m 2 - m' 2 ) 2 - Am 2 {k 2 - A 2 ) 

± k 



2(k 2 - A 2 ) 



2(k 2 - A 2 ) 



(4.20) 



There are three possible cases for Eq. ( |4.20| ): 



1. No solutions exist for Ex/2, or both solutions lie outside the integration interval 
[m, Ae\- Since there is at least one point (E = m), where the function in 
Eq. ([4.17 ) vanishes, we conclude that in this case the imaginary part disap- 
pears. 

2. One solution E\ lies inside the integration interval, the other one, E2 lies 
outside. It follows that the integral for the imaginary part has to be performed 
from Ex to A E in this case. 

3. The two solutions Ex < E 2 lie within the integration interval. In this case, the 
integral for the imaginary part has to be performed from Ex to E 2 - 

Since we will deal with integration limits of this type very often, we will adopt the 



convention that the solutions Ex/2 of Eq. ( |4.20| ) are tailored as described above. 
Eq. Q4.17 ) then takes the form 



dE f(±E - fx) log 



(A + E) 2 - (p - kf 



m 



(A + E) 2 - (p + k) 2 ~m 



a 



(4.21) 



+ i 7rSg ^ (A) [ E2 dE f(±E - d . 
k JEi 
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Together with Eq. flPOp , Eq. (|OlD provid es a sufficient basis for the numerical 
calculation of B^. It can be shown that in the limit k — > 0, Eq. (|4.21 ) continously ap- 



proaches Eq. (|4.13|) . Thus the generic function B^ is determined, and consequently 
B via Eq. Q. 

Figure |3] shows B for the same parameter set as was used in Fig. ^, but with 
k = lOOMeV. For k < k, Eq. (|4.20|) has one solution, from which one obtains an 
imaginary part. For ko > 2m, Eq. (|4.20| ) has two solutions, which become cut off at 
high ko, so that the imaginary part continously goes to zero. Figure ^ shows another 
special case. Here we plot B as a function of k for ko = 0. In this case, we have no 
imaginary part at all [|7|. 

D. Programming Considerations 

At this point, we deviate from the analytical calculation to discuss some of the 
basic aspects of programming such integrals. Comments made here are relevant for 
the B and Co integrals, the first of which is the least complicated integral containing 
multiple singularities. We discuss this now, bearing in mind that similar techniques 
are required for the more complex Co integral. 

To calculate Bo for the most general case via the decomposition ( |4.7| ), one has to 
take several things into consideration. First, the type of singularity determines the 
type of integration routine that should be utilized. One requires three such routines: 
(i) one routine which integrates smooth functions, (ii) one which integrates functions 
containing logarithmic or other integrable singularities, and (hi) one which computes 
the Cauchy principal value of an integral containing a 1/x singularity. Examples for 



this can be found in the literature 11.12 



On entry, the program has to decide, which of the forms (|4.13|) , pertinent to the 
case k = 0, or ( |4.21|) , pertinent for k > 0, to use. It is necessary to handle the 
two cases separately, since (i) Eq. ( [4.13| ) and Eq. ( f4.21| ) have a different analytical 



structure, and (ii) Eq. ( |4.21| ) contains a factor 1/k. In actual practice, the case 
k = is needed very often. After that, the position of singularities can be computed 
according to Eq. ( f4.11|) or Eq. (|4.20|) , respectively. The appropriate integration 



routine has to be chosen, depending on the result of this computation. 

Note that although the decomposition as given Eq. ( f4.7| ) is the most general 
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one possible, it does not always provide the optimal way of calculating B . If, for 
example, one knows a priori that an actual calculation will be confined to mi = m 2 , 
\i\ = ji2 and k = 0, the four terms of Eq. ( |4.7|) can be combined, yielding the 
simplified expression 

B (0, m,/i, m,fji,k , (3, A) = 8 v£ E dE(f(E -ft)- f(-E - //)) _ P ^ 2 (4.22) 

fc n \\ Jifcg - 4m 2 

9((2A E -A;o)(A ; o-2m)) . 

The result of this equation stays finite in the limit ko — > 0. The functions Bq, how- 
ever, individually diverge in this limit, so that a computation of B by decomposition 
suffers from cancellations at small ko- Note that similar simplifications follow for 
C . 




V. THE CALCULATION OF C 
A. Matsubara Summation 



In this section, we examine the analytical structure of the three fermion line 
integral Co, and decompose it into a sum involving a single generic function. 
The three fermion line function given in Eq. (|2.3|) has poles at 

iuo n = ± E\ (5.1a) 
iu n = iv m - fi 2 ± E 2 (5.1b) 
iio n = ia t - fx 3 ± E 3 (5.1c) 



with Ei, E 2 and E 3 as specified in Eq. ( |2.4| ). In analogy to Eq. ( |4.6| ), we define the 
complex variables 

Ai = iv m + Hi- H2 (5.2a) 
\ 2 = icti + fii - /j, 3 (5.2b) 
A 3 = iv m - i<*i ~ + = Ai - A 2 . (5.2c) 
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We again note that these parameters become real numbers after the analytical con- 
tinuation is performed. 

After summing the Matsubara frequencies, one obtains a sum over six terms that 
arise from each of the poles, 

Co(M,% ? , m 2 ,/i 2 ,w m , m 3 ,/j, 3 ,ia h (3, A) 



= 8tt 2 



d 3 p 
(2^)3 



E 1 l(\ 1 - El )i-El}i(\ 2 -E 1 y-Ei] 



+ (E 1 -> -^i) 



+ 



f(E 2 -» 2 ) 



E 2 [(Ax + £ 2 ) 2 - £ 2 ][(A 3 + E 2 f - E 2 3 ] 



+ (E 2 — > —E 2 ) 



+ 



f(E 3 - /i 3 ) 



E 3 [(A 2 + E 3 f - £ 2 ][(A 3 - E 3 f - El] 



(5.3) 



It is useful to make the substitution p — > A;— p in the third and fourth terms, p — ► <f— p 
in the fifth and sixth terms, in order to get rid of the angular dependence in the 
arguments of the Fermi functions. Again this is performed under the assumption 
that A ^> k, q. After this, we obtain a decomposition of Co of the form 

C (k,q,5^ m u n u m 2 ,/i 2 ,iz/ m , m 3 , /j, 3 ,ia h /3,A) 

= C^(-Ai, -A 2 , k,q,S^ TOi,m 2 ,m 3 ,/xi, (3, A) 
-Cq(\i, A 2 , k,q,S^ m 1 ,m 2 ,m 3 ,ii 1 , /3,A) 
+C , ( f(Ai,A 3 , fc, gl,^.^ m 2 ,mi,m 3 ,/i 2 , /3,A) 
-CS"(-Ai, -A 3 , k, \k- <z1,%_ 9 -, m 2 ,m 1 ,m 3 , /i 2 , /3,A) 
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+C h (A 2 , -A 3 , g, \q- k\,6fj_fr m 3 ,m 1 ,m 2 ,fi 3 , A A) 

-Co(-X 2 ,X 3 , q, \Q-%\,8f i f_fr m 3 ,mi,m 2 ,/i 3 , /3,A) , (5.4) 
where we have introduced the three fermion line generic function 

Cq(Xi,X 2 , k,q, 5, m, mi,m 2 , /i, /?, A) = 87r 2 (5.5) 
ci 3 p f(±E-n) 1 



x 



(27r) 3 £ [(Aj + £)2 _ (p - k)2 _ m 2] [ (A;j + £ )2 _ ^ 2 _ m 2] 



with = v^ 2 + m2 - in analogy to Eq. (4J3), one can derive a symmetry relation 
for the function Cq, as being 

C {k,q,5^ mi,/ii, m 2 ,fx 2 ,iu m , m 3 ,n 3 ,ia h /3,A) (5.6) 

= C (fc, | A; - g|,%^_^ ^2,^2, -iv m , m 3 , fi 3 ,iai - iv m , (3, A) . 

As was the case for the two line generic function Bq , C is calculated for several 
special cases. We list these in order of increasing difficulty: 

1. The case k = q = 0. 

2. The case k = 0, q > 0. 

— * 

3. The case k > 0, q > 0, p and A; collinear. 

4. The case k > 0, g > 0, p and & not collinear - general case. 

In the following sections, we will derive results for each of these special cases. 



B. Calculation of Cq for k = q = 



In this case, Eq. Q5.5p reduces to 

d 3 p f(±E-fi) 1 



Cf = 87T 2 



(2tt) 3 E [(X^ E) 2 -p 2 -m 2 ][(X 2 + E) 2 -p 2 -m 2 2 ] 



1 lim/ A ' d E J n±E -*l , (5.7) 



Ai A 2 6-^o 7 m (E-E x - ie) (E-E 2 - ie) 
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where we have again omitted the arguments of Cq for convenience and have included 
a small imaginary part in the denominators to do the analytical continuation. The 
constants E\ and E 2 are defined by 

E, = - A ' + ^" m ' (5.8a) 

If m < E\ = E 2 < Ae, the integral in Eq. (|5.7| ) diverges. We suppose therefore that 
Ei ^ E 2 , if one of them lies inside the integration interval. After taking the limit 
e — > 0, we obtain 

C± = A V f AE dE pf{±E-n) 

L [(\i+Ef-p*-mi][(\ 2 +E) 2 -p2-m$\ { ' 

+ l Xl x 2 {E 1 -E 2 ) ^ lf{±El ~ ^ ei{AE ~ El){El ~ m)) 

-p 2 f(±E 2 - /i)6((A i? - £ 2 )(£ 2 - m)) 

The O functions here guarantee that the imaginary part occurs only if the singu- 
larities appear inside the integration interval. As was detailed in Sec. |1V PL the 
numerical evaluation of the integral in Eq. ( |5.9| ) has to proceed differently if Ei 
or E 2 lie within the interval [m, Ae] (Cauchy integration), or if they lie outside of 
[m, Ae] (integration of a smooth function). 

Figure § shows m 2 Co for k = q = as a function of ko/A = 2qo/A for a similar 
parameter set as was used in Fig. 0. Again we obtain an imaginary part for ko > 2m. 



C. Calculation of Cq for k = 0, q > 

Without loss of generality, we assume that q points in the z direction. With this 
assumption, the scalar product pq occurring in Eq. ( |5.5| ) simplifies to pq cos 8. The 
integration over <p becomes trivial and one obtains 



Ct = - lim / d 3 p- 



x 



E Af + 2\iE + m 2 -m\- iesgn(Ai) 
1 

\ 2 + 2\ 2 E + 2pq cos 9 — q 2 + m 2 — m 2 — iesgn(A 2 
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1 



, lim 

2AiA 2 e-*0 



* dE pf{±E-n) /•+! 



dx- 



E-Ex-ieJ-i E - E 2 + f% - ie 



where now 



Eo 



A 2 , — q 2 + m 2 



2A 2 



(5.10) 



(5.11) 



and £?i is as defined in Eq. (|5.8aj) . The angular integral is singular if \X 2 {E — E 2 )\ < 
pq. This equation has exactly the same structure as Eq. ( f4.16| ), so Eq. ( |4.20| ) and the 
following discussion can be applied here too. We label the endpoints of the singular 
interval E 2 \ and E 22 . The angular integral becomes 



lim 



+ 1 



dx 



A, 



lim 



E - E 2 + - ie pq e->o J-i 



dx 



X 2 (E~E 2 ) 
pq 



— log 

pq 



\ 2 (E-E 2 ) 



I" I 



+ 1 



M{E-E 2 ) 
pq 



^, 7rA 2 sgn(A 2 ) Q / 



+ x - zesgn(A 2 ) 

X 2 {E — i? 2 ) 



pq 



log 



(A 2 + Ef - (p - qf 



m 



(A 2 + E) 2 — (p + q) 2 - m\ 
7rA 2 sgn(A 2 



+1- 



pq 



Q((E 22 - E)(E - E 2l )) 



(5.12) 



For we have accordingly 



1 r A E 
lim 



+ z 



2Aig e->o 
vrsgn(A 2 



log 



(A 2 + E) 2 -(p- qf 



(X 2 + Ef-(p + qf-m 2 



lim 

2Aig £^oJe 21 



(5.13) 



E - Ei - ie 

In this expression, we again have to take the limit e — * 0. The final expression we 
find is 



1 



2A 1(? 
7r 2 sgn(A 2 



E — Ei 



V dE 



log 



(A 2 + Ef -{p- qf 



(X 2 + Ef-(p + qf 



rrin 



+ i 



2Xiq 

nf(±Ei - y) 
2Xiq 



f(±Ei - m)©((£ 22 - Ei)(Ei - E 2 i)) 



log 



(A 2 + Eif - (pi - qf 



(A 2 + B1) 2 - (p, + qf - m| 



2Xiq Je 21 E- Ei 



e({A E -Ei)(Ei-m)) 

(5.14) 
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A numerical evaluation of Eq. ( 5.14 ) has to take into account Cauchy singularities 
at E — Ei and (integrable) logarithmic singularities at E — E 21 and E = E 22 . 



D. Calculation of Cq for k > 0, q > 0, q and k collinear 



Again, we assume without loss of generality that q and k have the form 







k 





\rjk J 



(5.15) 



where rj = ±1. With this assumption, we have pq = pqcosO and pk = rjpkcosd in 
Eq. (|5.5|). The integration is again trivial, and one obtains 



7T e^O I 



E Af + 2XiE + 2r]pk cos 6 — k 2 + m 2 — m 2 — iesgn(Ai) 
1 



A| + 2X 2 E + 2pq cos 6 — q 2 + m 2 — m 2 — iesgn(A2) 



-i-hm [ AE dEpf{±E-ri 



x dx 



1 E-E l + ^x-itE-E 2 + f 2 x-it 



(5.16) 



with 



Ex 



\{ — k 2 + m 2 — m{ 
2A^ 



(5.17) 



and E2 from Eq. ( |5.11| ). The integral over x becomes (cf. ||13|| , integral No. 12.8) 

1 1 



r+l 
lim / dx 



oy_! ~ E -E 1 + ^x-ieE-E 2 + f 2 X-ie 



A1A2 



j9gAi(E - £1) - ripk\ 2 (E - E 2 ) 

{\ 2 {E - £ 2 ) + pg)(A 1 ( J B - Ex) - V pk) 



x log 



(X 2 (E - E 2 ) - pq){\x{E - Ex) + rjpk) 
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+1 



7rAiA 2 r7sgn(Ai)0 (1 



Ai(g-gi) 
pk 



r]pkX 2 (E - E 2 ) -pq\i(E - E x 



7rAiA 2 sgn(A 2 )0 (l 



\ 2 {E-E 2 ) 
pq 



pqXx(E - Ex) - r]pk\ 2 (E - E 2 ) 
In what follows, we label the energies for which 



(5.1? 



\X 1 (E-E 1 )\ =pk 



(5.19) 



is fulfilled -Ei 1/2, the energies for which 



\X 2 (E - E 2 )\ = pq 



(5.20) 



is fulfilled E-ixii- These can be calculated explicitly via Eq. ( |4.20| ). Note that the 
discussion following Eq. (|4.20|) applies here too. 

The second and third terms of Eq. (|5.18|) still contain a singularity at E — Eq, 
where 



Eq = -(qXxEt - f]k\ 2 E 2 ) 



(5.21) 



with 



C = 2{qX l -T 1 kX 2 



(5.22) 



This is, however, not true for the first term, since the logarithmic factor goes to 
zero at E — Eq, making the integrand continous here. So we again add a small 
imaginary part in the denominators where appropriate and obtain 



C o = 7 / dE ~ p ^ lo § 

l, Jtn Hi — HiQ 



(X 2 (E-E 2 ) +pq){X 1 {E-E 1 )-r t pk) 



(X 2 (E - E 2 ) - pq)(Xi(E - Ex) + rjpk) 



.Tr^sgn(Ai) f(±E - n) , .Trs^Aa),. ,„ f{±E - /i) 



lim / dE- 

o J En E - Eq- ie 



+ i — - / — lim / dE- 

( e->o Je 2 i E — Eq — ie 



iE ^ ±E - ") io g 

Q Jm E — Eq 



(X 2 (E - E 2 ) + pq){Xx{E - Ex) - rjpk) 



(X 2 (E - E 2 ) - pq)(Xx(E - Ex) + rjpk) 



+ f(±E - /i)e((Ex2 - Eq)(E - Exx)) 
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7r 2 sgn(A 2 ; 
C 



■f(±E - fJ,)Q((E 22 - E )(E - E 21 )) 



7rr7sgn(Ai) r El2 , f(±E — a) 7rsgn(A 2 ) ^f E22 , f(±E — a) 

For the numerical evaluation of the real part, one has to take into account the 
logarithmic poles at E = E\\i 2 and E = E 2 \i 2 . For the imaginary part, the Cauchy 
singularity at E = E must be integrated. 

In Fig. H we show m 2 Co as a function of ko/A. for q = qo = ko/2 and all other 
parameters as in Fig. |5|. This special case has been considered previously in Ref. ||. 
As can be seen from Eq. (|5.4|), the computation of Co requires a computation of Cq 
for k = 0, q > 0, as well as for k, q collinear. 



E. Calculation of for the General Case 



In this subsection, we handle the general case of arbitrary momenta in the eval- 
uation of Cq from Eq. (|5.5| ). We choose the coordinate system in such a way that 





From this choice we obtain 



( ° ^ 

k sin 5 

i k cos 6 j 



(5.24) 



pk = pk(cos 9 cos 5 + sin 6 sin S cos </>) 



pq = pq cos 9 . 



It is useful to introduce the abbreviations 

\\ + 2Ai£ + m 2 - m\ - k 2 



Cl 



C-2 



2pk 

Xl + 2X 2 E + m 2 -m 2 -q 2 
2pq 



(5.25a) 
(5.25b) 

(5.26a) 
(5.26b) 



In terms of these variables, takes the form 
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4irkqe->oj m p Jo c 2 + cos# — 2esgn(A 2 J 

./o c\ + cos 6* cos 5 + sin # sin 5 cos — -ksgn(Ai) 

The additional difficulty introduced by the generality of handling arbitrary momenta 
is evident in the integration. As will be seen in the following, this can however 
still be dealt with analytically. 

To perform the <fi integration, we use the formula 

lim r = ; 271 (0(a 2 - 6 2 )sgn(a) + i®{b 2 - a 2 )) , (5.28) 

6^o7o a + bcos<p-ie \ a 2 _ 521 V / t> w , } 



which we prove in Appendix |A|. With this formula, one immediately obtains the 
result 

V±_ 1 [ Ae jvf(±E - fi) /•+ 1 ^e(A 1 )sgn(c 1 + xcos5)+«0(-A 1 )sgn(A 1 ; 



Of = — / dE^ ^ / dx 



2 kq Jm v •>-! + c 2 - «esgn(A 2 )) 

(5.29) 

where 

A x = ( Cl + xcos5) 2 -sin 2 5(l-a; 2 ) . (5.30) 

In the following subsections, we examine the remaining angular integrals that are 
required. 

1. The 9 Integration for Ai < 
In the next step, we calculate the integral 

r, x e(-*.)«°M . (5 . 31) 

J - 1 ^|A 1 |(x + c 2 -iesgn(A 2 )) 
We first note that for x — ±1, A x is always nonnegative. A x has zeros at 



x — xx — — Ci cos 5 — y 1 — c 2 sin 5 (5.32a) 

and 
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x = X2 = —c\ cos 5 + y 1 — c\ sin 5 (5.32b) 

with — 1 < a?i < X2 < +1. Consequently two conditions must be fulfilled in order to 
obtain Ai < 0: (i) |ci| < 1 and (ii) X\ < x < x 2 . From this it follows that we can 
apply the formula 

lim/ - = ^=(sgn(x + c)e(A)+*0(-A)) , 

\j {b - x){x - a){x + c - it) \J\&\ 

(5.33) 

where A = (a + c)(b + c), which is proven in Appendix [B|. Note that for A > 
sgn(x + c) is a constant in the interval [a, b], which is expressed by the symbol 
sgn(a; + c). 

Using Eq. (|5.33|) , one obtains 

-i 0(-A 1 )sgn(A 1 ) 



-i 



Ai| (a; + c 2 - zesgn(A 2 )) 



dx (5.34) 



0(1 _ c 2 f s Sn(Xi) (ggn(c2 _ ^ cos5)Q(Ao) + , sgn(A2 )0 ( _ Ao )) 
|A | 



where the factor sgn(x + c 2 ) has been taken at x = (x\ + x 2 )/2 = —c\ cos 5. Ao can 
be obtained from Ai by substituting — c 2 for x: 

Aq = c\ + c 2 — 2cic 2 cos 5 — sin 2 5 (5.35) 



2. The 9 Integration for Ai > 

It remains to compute the integral 
f +1 6(Ai)sgn(ci + xcos5) , ^f +1 @(Ai)sgn(ci + x cosS) , 

llm / nr- / : TTT^dx = V 7=7 r dx 

e-*oJ_i ^/A^x + c 2 - zesgn(A 2 J) -/-l vAi(z + c 2 ) 

+ / sgn(A 2 )sgn(c i -c 2 cos^) e(Ao)e(i _ c , ) (g ^ 
v Ao 

In doing so, one easily notes that C\ +xcos<5 vanishes only at points where Ai < 0, 
so that the sgn term can be treated effectively as a constant. The real part can 
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be computed from [14], integral No. 231.10. Bearing the definition of X\i 2 from 
Eq. ( |5.32|) in mind, one obtains 

0(Ai)sgn(ci + xcos<5) 
VI 7=-, ; dx 



'A 1 (x + c 2 

e(A 



sgn(ci + cos S) 



F(+l) - 0(1 - c\)F{x 2 ) 



— sgn(ci — cos S) 



F(-l)-Q(l-c 2 1 )F(x 1 



e(-A ) 

' /=A 



L 



sgn(ci + cos S) 



G(+1)-G{x 2 ) 



— sgn(ci — cos S) 



G(-l) - G(x\) 



where 



G(x) = arccos 



(c 2 - ci cos 5) (x + c 2 ) - A c 



sin 5y 1 — cf \x + c 2 | 

and F is either one of the functions 

(x + c 2 )(ci cos 5 - c 2 ) + A =F VAqAi 



F+(x) = ±1oe 



X + c 2 



Eq. ( |5.37|) can be further simplified by 



F±(x 1 ) = F±(x 2 ) = ±log sin 5Jl - cf 



G(xi) = G(x 2 ) = 



(5.37) 



(5.38) 



(5.39) 



(5.40a) 



(5.40b) 



Up until this point, the calculation in this section has been purely analytical. 
However, although F + and F_ are equivalent in specifying the indefinite integral 
pertinent to Eq. ( |5.36| ), they are not equally well behaved numerically if either |ci| 
or |c 2 | is approximately equal to one. It is therefore necessary at this point to devise 
a numerically stable procedure, and we give our algorithm here. To this end, it is 
useful to define H as 



sgn(ci + cos 5) 



— sgn(ci — cos 5) 



F(+l)-Q(l-cl)F(x 2 ) 



F(-l)-e(l-4)F( Xl ) 



(5.41) 
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Using the identity 

F+(x) -F-(x) = log (sin 2 5 |l - c?|) , 
one can show that the following algorithm gives the correct result: 

• If | |c 2 | — 1| < | |ci| — 1|: Compute S as 

S = sgn(ci + cos5)F„(+l) — sgn(ci — cos5)F_(— 1) 

+ 6(cos 2 5 - c 2 )sgn(cos5) log sin 2 5(1 - c 2 ) 

• If ||c 2 | - 1| > ||ci| - 1|: 

— If c 2 — Ci cos 5 > 0: Compute S as 

S = sgn(cx + cos5)F + (+l) — sgn(ci — cos5)F_(— 1) 

— 0(c 2 — cos 2 5)sgn(ci) log sin 2 5(1 — c 2 ) 

— If c 2 — Ci cos 5 < 0: Compute S as 

S = sgn(ci + cos5)F_(+l) — sgn(cx — cos5)F + (— 1) 

+ 0(c 2 — cos 2 5)sgn(ci) log sin 2 5(1 — c\) 
This is implemented in the numerical procedure. 

3. Putting the Parts together 
The final expression for Cq is now given. We may write 



1 /(±S_^) 



f 'rfE JV ^ (e(Ao)S-e(-Ao)7rsgn(A 1 A 2 ) 

+6(-A )(sgn(ci + cos5)G(+l) - sgn(ci - cos<J)G(-l) 
+i7r9(A o )(sgn(Ai)0(l - c 2 )sgn(c 2 - c x cos 5) 

+sgn(A 2 )0(l - c 2 )sgn(ci - c 2 cos 5)) 
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where H has to be taken from Eq. Q5.43 ). In the spirit of Sec. IV D , in order to 



evaluate this integral, one has to identify the singular points of the integrand. It is 
easy to see that the following points are singularities: 

• Energies, at which |ci| = 1 or | C2 1 = 1. These can be calculated from Eq. fl4.2U| ). 

• Energies, at which Ao = 0. To compute these energies, one notes that c\ and 
c 2 can be written as 

ai + b\E a 2 + b 2 E 

ci = c 2 = , (5.45) 

p p 

where a, and 6j can be obtained from Eq. ( |5.26| ). With these constants, the 
equation A = can be cast into the form 

= E 2 {b\ + b\- 26x62 cos 5 - sin 2 6) 
+ 2E(aibi + a 2 b 2 - {a\b 2 + a 2 &i) cos 5) 
+ (a 2 + a 2 — 2a x a 2 cos 5 + m 2 sin 2 5) (5.46) 

which can be easily solved for E. 

• Energies, at which |ci| = | cos<5|. These can be computed by 

E= ^m h 72V (^ a i ± k cos Sy/al - 4m 2 (A 2 - k 2 cos 2 5)) . (5.47) 

Z[rC COS — ^lj ^ ' 

All of these points correspond to integrable singularities. On closer inspection, one 
also obtains that the energies, at which the arguments of the logarithmic function (in 
F±) vanish, or the energies, at which the argument of the arccosine (in G) becomes 
greater than one, are included in these cases. 

To illustrate C for the general case, we plot m 2 C in Fig. [F|as a function of k /A 
with k = 3/4fc , q = 3/4g and S = 71/6. All other parameters are as in Fig. ||. 

VI. APPLICATION: KAON MASSES IN THE NJL MODEL 

In this section, we give a specific example for a calculation using our elementary 
integrals. It is also intended for the benefit of our reader, who wishes to compare 
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his/her results with our numbers. To this end, we compute the kaon masses in 
the framework of the SU(3) Nambu-Jona-Lasinio model ||,f||| as a function of 
temperature. In this model, the quark masses have to be determined from the 
coupled gap equations 

m q = m 0q + AGNjX^S^x, x) + 2KN 2 c (iti 1 S q (x, x))(itr 7 S s (a;, x)) (6.1a) 

m s = m0s +AGNJ,tx 7 S s {x,x) + 2KN 2 (itr 1 S q (x,x))(itr 1 S q (x,x)) (6.1b) 

with the current quark masses mo q , m Qs , the number of colors N c and the coupling 
constants G and K. S* denotes the finite temperature quark propagator 

Sf{x-x\r-T')= l -Ye-^-^ f^-— r ^ x ~ x) (62) 

The trace of the propagator can be expressed in terms of the function A: 

itr 7 S^(x, x) = —^^A(mf,jj,f) (6.3) 
and the gap equation takes the form 

m q = m 0q - ^m q A(m q , fj, q , (3, A) (g - ^^m s A(m s , /i s , /3, A)^ (6.4a) 

GN KN 2 
m s = m 0s -^-m s A(m s , /i s , (3, A) + — (m q A(m q , /j, q , ft, A)) 2 (6.4b) 

Meson masses are computed from the quark-antiquark scattering matrix. The 
scattering matrix for light and strange quarks in the pseudoscalar channel can be 
written as 

M = i-J<% ,k) • (6 - 5) 

where (h®, k) is the total four momentum of the quarks, Tl(ko, k) the irreducible 
polarization and K e R an effective coupling constant, which is given by 



KN C KN C 
K cS = G H j^itr^S q (x, x) = G - -^-m q A(m q , fi q , (3, A) (6.6) 



The irreducible polarization II(A;q, k) is computed from the diagram in Fig. [8|. 
One obtains 
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-ill(iv m , k) = -n c - 



d 3 p 



tr- 



iS q (iuj n ,p)ij 5 iS s (iuj n - iv m ,p- k)i^ b 



I ( iuJn + Li<q)( iuJ n - iv m + fx s ) -pjp-k) -m q m s 

1 C P^J (27T)3 W iUn + jIq y_ E 2\^ iuJn _ iUm + fMs) 2_ E 2 ] 



(6.7) 



where £7 g = yp 2 + m 2 and E s = y(p—k) 2 + m? s . To cast this into a form which 
contains the functions A and Bo, one uses the identity 



(iu n + fi q )(iu n - iv m + fjL s ) - p{p- k) - m q m s 



iu n + fi q ) 2 - E 2 ) + ((iu; n - %v m + Us) 2 - 



- (fig -fi s + w m ) 2 + k' 
After continuing iv m to k , one immediately obtains the result 



+ (m a - m s 



(6. 



n(fc ,fc) 



A(m q ,(i q , p,A)+A(m a ,ii s , (3, A) 



(6.9) 



+ ((m, - m s ) 2 - (fco + /ig - /x s ) 2 + A; 2 ) B (k, m q , fj, q , m s , /i s , k , (3, A) 
The kaon mass is now computed using the dispersion relation 



1 - 4K eS Uim K , 0) = 



(6.10) 



Using the numbers m 0q = 5.5MeV, m 0s = 140.7MeV, fj, q = = 0, A = 602.3MeV, 
GA 2 = 1.835 and KA 5 = 12.36, one obtains the zero temperature results m q = 



367.7MeV and m s = 549.5MeV from Eq. ( |6.4| ). Solving Eq. ( |6.10|) with these num- 
bers yields m K = 497.7MeV. The temperature dependence of m q , m s and m K is 
computed using the same formulae and the result is depicted in Fig. []. Further 
examples can be found in Ref . |§ . 



VII. SUMMARY AND CONCLUSION 



In this paper, we have presented the technical aspects required for a numerical 
evaluation of the one, two and three fermion line integrals at finite temperature, 
structuring each as a sum of terms of a specific generic integral taken at different 
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values of its arguments, that are often required for field theoretic calculations such 
as within the Nambu-Jona-Lasinio model. Both real and imaginary parts of these 
functions are explicitly calculated. Concomitantly, we have illustrated all functions 
graphically for certain parameter sets, in order to enable the user to verify his/her 
own calculation. A simple example, employing two of the integrals has also been 
shown. A computer source code is also available that routinizes these integrals. The 
generalization to bosonic line integrals can also be seen to follow in an analogous 
fashion. 
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APPENDIX A: PROOF OF EQ. ( gjgg ) 



Our assertion is 

limT , =-. 271 (e(a 2 -b 2 )sgn(a)+z&(b 2 -a 2 )) . (Al) 



<^o Jo a + b cos <p — ie yja^ — b 2 
The real part of this equation can be easily obtained from integral tables. (See e. g 



141 ) integral No. 331.41.) The singularity in the integral only appears for \b\ > \a\, 
which gives the function for the imaginary part. In this case, singularities appear 
at 

= arccos(— a/b) and <fi = 2tt — arccos(— a/b) . (A2) 

The first of these singularities appears in the interval [0, 7r], the second one in the 

interval [tt, 2tt] . We split the integral into two parts 

^ d<j> 1 
hm / = hm — 



e^o Jo a + b cos <p — ie e->o b Jo a/b + cos <p — iesgn(6) 



, 1 d( t> /A N 

+ lim T/ -71 T— 777 A3 

e^o b J tt a/b + cos0 — tesgn{b) 
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and focus on the first of these two parts. Let O = arccos(— a/b) and the function 
/(0) be defined by 

{(cos + a/b) /((f) - O ) for 7^ 0o 
(A4) 
sin 0o = y^I — (a/b) 2 for = O 

It can be easily seen that /(0) is continous in the interval [0, tc] and 

1 Z" 71 " dcj) . 1 r 

lim - / — ■ — = hm ■ 



o b Jo a/b + cos0 — iesgn(b) b Jo (0 — 0o)/(0) — iesgn(b) 

lim - r d(t> (A5) 

e-o&io (0-0 o -^sgn(6))/(0) 1 ] 



since /(0o) is positive. Our standard formula Eq. ( (4.12|) can be applied to this 
integral to yield 

Ve-ofeio a/6 + cos0-zesgn(6); b & V ; /(0 O ) - a 2 

The second integral on the right hand side of Eq. (|S3|) can be treated in the same 
fashion and is found to give the same contribution. QED. 

APPENDIX B: PROOF OF EQ. 

The assertion is 

lim [ b j= . . = — = -jL= (sgn(x + c)0(A) + iQ(-A)) (Bl) 



' a — x) (x — ~a)(x + c — ie) y\A 
with A = (a + c){b + c). To prove this, we first apply Eq. (|4.12|) to obtain 



f b dx r b dx in<d(—A) 
lim / = = VI = + 



e-*0J a 



^(b — x)(x — a)(x + c — ie) ^° JJb — x)(x — a)(x + c) V~ A 

(B2) 



The G function on the right hand side of Eq. ( |B2| ) emerges from the fact that the 
singularity occurs only if 

a <- c <b <=> A < . (B3) 
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The real part of Eq. (|B2|) can be looked up in integral tables. (See e. g. fll3| , integral 



No. 221.7.) For A > 0, one obtains 



dx 1 ( — (c + a)(b — x) + (c + b)(x — a)\ , 
arcsin v A \7 \ " ( B4 ) 



yj(b-x)(x-a)(x + c) V (6-a)|x + c 

which gives 



dx 7rsgn(x + c) ^ B5 ^ 



' a ^{b-x){x- a){x + c) 
For A < 0, one obtains 

dx 1 , (J (c + a)(b-x)- J-(c + b)(x- a)Y 

log 



yj(b - x){x - a)(x + c) v 7 -^" 



\x + c 



(B6) 

which gives 

Vf b dX = . (B7) 

V ~ x ) ( x — a ) ( x + c ) 

Combining Eqs. (P2]), flBli| ) and QB"7|), proves the assertion. 
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FIGURES 



FIG. 1. A/ A 2 as a function of m/A at fi = 0, (3 -► oo and A = 602.3MeV. 

FIG. 2. S for fc = 0, m = fi 2 = 0, m x = m 2 = 367.7MeV, T = and A = 602.3MeV. 
The solid line gives the real part, the dashed line the imaginary part. 

FIG. 3. B for k = lOOMeV, m = fi 2 = 0, mi = m 2 = 367.7MeV, T = and 
A = 602.3MeV. The solid line gives the real part, the dashed line the imaginary part. 

FIG. 4. B for k = 0, m = fi 2 = 0, m 1 =m 2 = 367.7MeV, T = and A = 602.3MeV. 
Note that Bq is a real function in this case. 

FIG. 5. m 2 C for k = q = 0, mi = m 2 = m 3 = 367.7MeV, m = fi 2 = fi 3 = 0, 
qo = 2ko, T = and A = 602.3MeV. The solid line gives the real part, the dashed line 
the imaginary part. 

FIG. 6. m 2 Co for k = 0, q = qo = 2fco, m\ = m 2 = 777,3 = 367.7MeV, n\ = fi 2 = /J3 = 0, 
T = and A = 602.3MeV. The solid line gives the real part, the dashed line the imaginary 
part. 

FIG. 7. 77i 2 C for k = S/4k , q = 3/4q , 5 = vr/6, m 1 = m 2 = m 3 = 367.7MeV, 
Mi = A*2 = A*3 = 0, T = and A = 602.3MeV. The solid line gives the real part, the 
dashed line the imaginary part. 

FIG. 8. Feynman diagram for the irreducible polarization function. 

FIG. 9. Temperature dependence of 777 g (solid line), tt7 s (dashed line) and mx 
(dot-dashed line). 
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